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Summary. The discovery of an excess of light in the far-ultraviolet (UV) spectrum 
in elliptical galaxies was a major surprise in 1969. While it is now clear that this 
UV excess is caused by an old population of hot helium-burning stars without large 
hydrogen-rich envelopes rather than young stars, their origin has remained a mys- 
tery. Here we show that these stars most likely lost their envelopes because of binary 
interactions, similar to the hot subdwarf population in our own Galaxy. This has 
major implications for understanding the evolution of the UV excess and of ellipti- 
cal galaxies in general. In particular, it implies that the UV excess is not a sign of 
age, as had been postulated previously, and predicts that it should not be strongly 
dependent on the metallicity of the population. 



1 Introduction 

One of the first major discoveries soon after the advent of UV astronomy was 
the discovery of an excess of light in the far-ultraviolet (far-UV) in elliptical 
galaxies (see the review by O'Connell PQ). This came as a complete surprise 
since elliptical galaxies were supposed to be entirely composed of old, red 
stars and not to contain any young stars that radiate in the UV. Since then 
it has become clear that the far-UV excess (or upturn) is not a sign of active 
contemporary star formation, but is caused by an older population of helium- 
burning stars or their descendants with a characteristic surface temperature 
of 25,000 K [2], also known as hot subdwarfs. 

The origin of this population of hot, blue stars in an otherwise red pop- 
ulation has, however, remained a major mystery [3J. Two scenarios, referred 
to as the high- and the low-metallicity scenario, have been advanced. In the 
low- metallicity model [3] , it is argued that these hot subdwarfs originate from 
a low-metallicity population of stars which produce very blue helium core- 
burning stars. This model tends to require a very large age of the population 
(in fact, larger than the generally accepted age of the Universe); it is also not 
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clear whether the population is sufficiently blue to account for the observed 
UV color. Moreover, the required low metallicity appears to be inconsistent 
with the large metallicity inferred for the majority of stars in elliptical galax- 
ies [5]. In contrast, the high- metallicity model [5] [7] assumes a relatively high 
metallicity - consistent with the metallicity of typical elliptical galaxies (~ 1 - 

3 times the solar metallicity) - and an associated enhancement in the helium 
abundance and, most importantly, postulates an enhanced and variable mass- 
loss rate on the red-giant branch, so that a fraction of stars lose most of their 
hydrogen-rich envelopes before igniting helium in the core pj [5] . 

Both models are quite ad hoc: there is neither observational evidence for 
a very old, low-metallicity sub-population in elliptical galaxies, nor is there 
a physical explanation for the very high mass loss required for just a small 
subset of stars. Furthermore, both models require a large age for the hot com- 
ponent and therefore predict that the UV excess declines rapidly with redshift. 
This is not consistent with recent observations, e.g. with the Hubble Space 
Telescope (HST) [9]. In particular, the recent survey with the GALEX satel- 
lite [10] showed that the UV excess, if anything, may increase with redshift. 
Indeed, the wealth of observational data obtained with GALEX is likely to 
revolutionize our understanding of elliptical galaxies. While Burstein et al. 
[11] appeared to have found a correlation between the UV-upturn and metal- 
licity in their sample of 24 quiescent elliptical galaxies, which could support 
the high-metallicity scenario, this correlation has not been confirmed in the 
much larger GALEX sample [10], casting serious doubt on this scenario. 

Both models ignore the effects of binary evolution. On the other hand, hot 
subdwarfs have long been studied in our own Galaxy [T2|, and it is now well 
established [13] that the vast majority of (and quite possibly all) Galactic hot 
subdwarfs are the results of binary interactions where a star loses all of its 
envelope near the tip of the red-giant branch by mass transfer to a companion 
star or by ejecting it in a common-envelope phase, or where two helium white 
dwarfs merge to produce a single object (see [HUH] for references and details). 
In all of these cases, the remnant star ignites helium and becomes a hot 
subdwarf. The key feature of these binary channels is that they provide the 
missing physical mechanism for ejecting the envelope and for producing a hot 
subdwarf. Moreover, since it is known that these hot subdwarfs provide an 
important source of far-UV light in our own Galaxy, it is not only reasonable 
to assume that they will also contribute significantly to the far-UV in elliptical 
galaxies, but is in fact expected. 

2 The Model 

To quantify the importance of the effects of binary interactions on the spectral 
appearance of elliptical galaxies, we have performed the first population syn- 
thesis study of galaxies that includes binary evolution (see also [TBI [T71 ITS] ). 
It is based on a binary population model [HI [To] that has been calibrated 
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to reproduce the short-period hot subdwarf binaries in our own Galaxy that 
make up the majority of Galactic hot subdwarfs [T3]- The population synthe- 
sis model follows the detailed time evolution of both single and binary stars, 
including all binary interactions, and is capable of simulating galaxies of arbi- 
trary complexity, provided the star-formation history is specified. To obtain 
galaxy colors and spectra, we have calculated detailed grids of spectra for hot 
subdwarfs using the ATLAS9 [19] stellar atmosphere code. For the spectra and 
colors of single stars with hydrogen-rich envelopes, we use the comprehensive 
BaSeL library of theoretical stellar spectra [2D1 [5T] . 

3 Results and Discussion 




Fig. 1. The evolution of the far-UV spectrum with time for a single population 
where all stars formed at the same time. The flux fx is scaled relative to the visual 
flux (/ v ). 

Figure [1] shows our simulated evolution of the far-UV spectrum of a galaxy 
in which all stars formed at the same time, where the flux has been scaled 
relative to the visual flux (between 5000 and 6000A) to reduce the dynam- 
ical range. At early times the far-UV flux is dominated by the contribution 
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Fig. 2. Evolution of far-UV properties [the slope of the far-UV spectrum, /3puv, 
versus (1550 — V)] for a two-population model of elliptical galaxies. The age of 
the old population is assumed to be 12 Gyr (filled squares, filled triangles, or filled 
circles) or 5 Gyr (open squares, open triangles, or open circles). The mass fraction 
of the younger population is denoted as / and the time since the formation as t m inor 
[plotted in steps of Alog(t) = 0.025]. Note that the model for / = 100% (the dotted 
curve) shows the evolution of a simple stellar population with age i m inor. The legend 
is for &fuv, which is the fraction of the UV flux that originates from hot subdwarfs 
resulting from binary interactions. The effect of internal extinction is indicated in the 
top-left corner, based on the Calzetti internal extinction model with E(B — V) — 0.1 
[27] . For comparison, we also plot galaxies with error bars from HUT [22] and IUE 
observations [11] , The galaxies with strong signs of recent star formation are denoted 
with an asterisk (NGC 205, NGC 4742, NGC 5102). 



from single young stars. Binary hot subdwarfs become important after about 
1.1 Gyr, which corresponds to the evolutionary timescale of a 2M star and 
soon start to dominate completely. After a few Gyr the spectrum no longer 
changes appreciably. 

There is increasing evidence that many elliptical galaxies had some re- 
cent minor star-formation events [UJ HI] i which also contribute to the far-UV 
excess. To model such secondary minor starbursts, we have constructed two- 
population galaxy models, consisting of one old, dominant population with an 
assumed age t id and a younger population of variable age, making up a frac- 
tion / of the stellar mass of the system. In order to illustrate the appearance 
of the galaxies for different lookback times (redshifts) , we adopted two values 
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for toid, °f 12 Gyr and 5 Gyr, respectively; these values correspond to the ages 
of elliptical galaxies at a redshift of and 0.9, respectively, assuming that the 
initial starburst occurred at a redshift of 5 and adopting a standard ylCDM 
cosmology with H = 72km/s/Mpc, J?m = 0.3 and Qa — 0.7. Our spectral 
modelling shows that a recent minor starburst mostly affects the slope in the 
far-UV SED. We therefore define a far-UV slope index /3 F uv as fx ~ A* uv , 
where /3fuv is fitted between 1075A and 1750A. This parameter was obtained 
from our theoretical models by fitting the far-UV SEDs and was derived in 
a similar manner from observed far-UV SEDs of elliptical galaxies [TTJ [22] , 
where we excluded the spectral region between 1175A and 1250A, the region 
containing the strong Lya line. In order to assess the importance of binary 
interactions, we also defined a binary contribution factor b — -Fb/-Ptotal, where 
Fb is the integrated flux between 900A and 1800A radiated by hot subd- 
warfs produced by binary interactions, and F to tai is the total integrated flux 
between 900A and 1800A. Figure shows the far-UV slope as a function of 
UV excess, a potentially powerful diagnostic diagram which illustrates how 
the UV properties of elliptical galaxies evolve with time in a dominant old 
population with a young minor sub-population. For comparison, we also plot 
observed elliptical galaxies from various sources. Overall, the model covers 
the observed range of properties reasonably well. Note in particular that the 
majority of galaxies lie in the part of the diagram where the UV contribution 
from binaries is expected to dominate (i.e. where b > 0.5). 

The two-component models presented here are still quite simple and do 
not take into account, e.g., more complex star-formation histories, possible 
contributions to the UV from AGN activity, non-solar metallicity or a range 
of metallicities. Moreover, the binary population synthesis is sensitive to un- 
certainties in the binary modelling itself, in particular the mass-ratio distri- 
bution and the condition for stable and unstable mass transfer [T5]. We have 
varied these parameters and found that these uncertainties do not change the 
qualitative picture, but affect some of the quantitative estimates. 

Despite its simplicity, our model can successfully reproduce most of the 
properties of elliptical galaxies with a UV excess: e.g., the range of observed 
UV excesses, both in (1550 — V) and (2000 — V) (e.g. [25]), and their evolution 
with redshift. The model predicts that the UV excess is not a strong function 
of age, and hence is not a good indicator for the age of the dominant old 
population, as has been argued previously [55], but is very consistent with 
recent GALEX findings [10]. We typically find that the (1550 - V) color 
changes rapidly over the first 1 Gyr and only varies slowly thereafter. This 
also implies that all old galaxies should show a UV excess at some level. 
Moreover, we expect that the model is not very sensitive to the metallicity of 
the population since metallicity does not play a significant role in the envelope 
ejection process (although it may affect the properties of the binary population 
in more subtle ways). 

Our model is sensitive to both low levels and high levels of star formation. 
It suggests that elliptical galaxies with the largest UV excess had some star 
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formation activity in the relatively recent past (~ 1 Gyr ago). AGN and su- 
pernova activity may provide supporting evidence for this picture, since the 
former often appears to be accompanied by active star formation, while su- 
pernovae, both core collapse and thermonuclear, tend to occur mainly within 
1-2 Gyr after a starburst in the most favoured supernova models. 

The modelling of the UV excess presented in this study is only a starting 
point: with refinements in the spectral modelling, including metallicity effects, 
and more detailed modelling of the global evolution of the stellar population 
in elliptical galaxies, we suspect that this may become a powerful new tool 
helping to unravel the complex histories of elliptical galaxies that a long time 
ago looked so simple and straightforward. 

This work has been supported by the Chinese National Science Foundation 
under Grant Nos. 10433030 and 10521001 (ZH). See [28] for a detailed version 
of the paper. 
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